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ABSTRACT 

We present global structural parameter measurements of 109,533 unique, i?Fi60W-selected objects 
from the CANDELS multi-cycle treasury program. Sersic model fits for these objects are produced 
with GALFIT in all available near-infrared filters (-HF160W, ^Fi25W an d, for a subset, F F10 5w)- The 
parameters of the best-fitting Sersic models (total magnitude, half-light radius, Sersic index, axis ratio, 
and position angle) are made public, along with newly constructed point spread functions for each 
field and filter. Random uncertainties in the measured parameters are estimated for each individual 
object based on a comparison between multiple, independent measurements of the same set of objects. 
To quantify systematic uncertainties we create a mosaic with simulated galaxy images with a realistic 
distribution of input parameters and then process and analyze the mosaic in an identical manner as 
the real data. We find that accurate and precise measurements - to 10% or better - of all structural 
parameters can typically be obtained for galaxies with 7?fi60W < 23, with comparable fidelity for 
basic size and shape measurements for galaxies to -ffFi60W ~ 24.5. 



1. INTRODUCTION 

The Cosmic Assembly Near-Infrared Deep Extragalac- 
tic Legacy Survey (CANDELS) constitutes 902 or- 
bits of Hubble Space Telescope (HST) observing time. 
One of the primary motivations for CANDELS is the 
investigation of galaxies at redshifts z ~ 1.5 — 3, 
and in par ticular their struc tural and morphological 
properties (jGrogin et al.l 120 111) . This paper describes 
the characterization of structural gal axy properties in 
the H ST WFC3/IR imaging mosaics (|Koekemoer et al.l 
1201 lh through fitting the observed surface brightness 
distributions by two-dimensional parametrized models, 
whose surfac e brightness profiles follow the Sersic law 
(|SersidlT968l ). The resulting catalogs are made public 
online. 

Over the previous decade WFPC2 and ACS have en- 
abled comprehensive structural and morphological stud- 
ies of the galaxy population up to z ~ 1 at 'opti- 
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cal' wavelengths (e.g.. iGiavalisco et al.l [2004; Ri x et al.1 
12001 IScoville et aljl2007h . At redshifts z > 1 WFPC2 
and ACS observations sample the rest-frame ultra-violet, 
which complicates the interpretation of galaxy images, 
as dust can strongly attenuate the light and young stars 
that contribute little to the underlying stellar mass dis- 
tribution dominate over the older population. 

For about a decade now, deep near-infrared sur- 
veys have been conducted to remedy this (e.g., 
Labbe et all 120031: iFontana et al.l 120041 iQuadri et al.1 
20071: lLawrence et al.l 120071; iCapak elaLl I2007D . These 
ground-based surveys have provided large samples, and 
despite the limited spatial resolution (typically 0.4-0.8", 
or ~ 3 — 6 kpc at z > 1), general and fundamen- 
tal struc tural properties have bee n measured out to 
z ~ 2 (iTruiillo et all l200l 120061; iFranx et al.l l200l 
I Williams et ah 2009; C hang et al.ll2012T ). Galaxies were 
shown to be smaller in the past, and at all redshifts their 
sizes were found to correlate with star-formation activity. 
As in the local universe, galaxies with low star-formation 
activity tend to be smaller than galaxies with high (or 
normal) star-formation rates. 

The use of ground-based near-infrared imaging for the 
purpose of investigating the internal structure of galaxies 
and its evolution has largely been limited to simple size 
measurements. Any examination of galaxy structure be- 
yond this requires HST resolution. Near-infrared obser- 
vations with NICMOS over relatively small areas and tar- 
geted sampling of small numbers of pre-selected galaxies 
have before revealed some fundamental aspects of struc- 
tural properties at z ~ 2. Beyond confirming the size 
evolution of galaxies, the profoundly different structure 
of high-z galaxies could now be fully apprec iated. The 
discovery of mas sive, yet very small, galaxies (jZirm||2007l ; 
iToft et aT1l2007l ) posed a surprise to the community, and 
has instigated much debate regarding the formation and 
evolution of the most massive galaxies. Moreover, the 
concentration of z ~ 2 galaxy light profiles, presum- 
ably tracing the bulge-to-disk ratio, was found to corre- 
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late with star-formation activity, indicating that galaxies 
with low s tar-formation rate s have larger bulges at all 
z < 2 fe.g.. lKrieket~aT1[2009l) . 

The arrival of WFC3 has now allowed the exploration 
of galaxy structure at z > 1 with unprecedented data 
quality a nd sample sizes. The Ultra Deep Field (UDF) 
program ( Bouwcns et al.l I2010D broke ground in terms 
of depth and resolution, confirming previous claims on 
the structural properties of z ~ 2 galaxies and revealing 
further detail on the ir morphology and structure (e.g., 
iSzomoru et all l2010| ). Results from the HST/WFC3 
Early Release Science (ERS) program ((Windhors t et al.l 
1201 If) foreshadowed the power of CANDELS. That mod- 
estly sized ERS program provided data quality and quan- 
tity on par with the largest NICMOS data sets i n ex- 
istence (|Scoville etHI 120071 : iConselice et al.l 1201 If ) and 
produced new insights into the ba s ic structural param - 
eters (e.g.. Ivan der Wei et al.ll2011t iCassata et al.ll20Tll ) 
and morphologies fe.g.. iCameron et al. 1 1201 ID of z ~ 2 
galaxies. 

CANDELS fills the gap between the ground-based sur- 
veys and the existing HST near-infrared surveys. It pro- 
vides much larger samples than previous near-infrared 
HST data sets and much improved depth and res o lution 
co mpared to ground -based surveys. iWuvts et al.1 (|201 lh 
and lBell et alj (|2012f ) showed with unprecedented clarity 
how star-formation activi ty an d structure are st rongly 
related. iPapovich et~aT1 (|2012f ) and lLotz et al.l (|201lf ) 
examined the environmental dependence of galaxy sizes 
and merger frequency, respectively. Moreover, the larger 
area allows to study the p roperties of more rare ob- 
jects. iKocevski et a l. (2012) compared the morphologies 
of galaxies that do and do not host AG N and concluded 
there is no significant difference, and iKartaltepe et al.l 
(|2011f) found a tentative connection between extreme star 
formation activity and merging. 

CANDELS also allows us to probe galaxy structure 
down to previously unattainable low mass and luminos- 
ity limits - apart from the UDF - in particular through 
the deep segment of the survey. Furthermore, the struc- 
ture of massive galaxies beyond z = 3, where knowledge 
is still sparse, c an be explored (e.g.. ICaputi et al1l2012t 
IGuo et al.|[2012f ). Finally, CANDELS aims at obtaining 
for the first time a comprehensive view of the morpholo- 
gies of z > 4 Lyman-break and Ly-a emitters at wave- 
lengths longer than 2000A in the rest frame. 

In this paper we describe the measurements of struc- 
tural parameters of 109,533 unique objects in the CAN- 
DELS WFC3/IR data, representing roughly 2/3 of the 
full survey. Our online materials will be updated with 
the final 1/3 of the survey once observations have been 
completed by the end of 2013. We describe the imag- 
ing, object detection, and background flux level esti- 
mation in Jj2] In <j3] we construct Point Spread Func- 
tion (PSF) models that are used as the default PSFs 
in the CANDELS collaboration; t hese models a r e also 
made public. Using GALAPAGOS (jBarden et al.l l2012f ) 
we prepare flux and noise im age cutouts and describe 
how GALFIT (|Peng et al.ll2010|) is used to produce single- 
component Sersic model fits (@. We estimate random 
and systematic uncertainties in the parameter estimates 
through the comparison between different image data 
sets of the same objects, and Sersic fits to simulated 



Table 1 

Publication dates of CANDELS catalogs with GALFIT -based 
structural parameters published online. Those in bold face are 
published along with this paper. The rest will be published over 
the next year, within 4 months after the last observation of each 

field. 



Field 


F105W 


F125W 


F160W 


COSMOS 




12/2012 


12/2012 


EGS 




09/2013 


09/2013 


GOODS-N 


12/2013 


12/2013 


12/2013 


GOODS-S 


12/2012 


12/2012 


12/2012 


UDS 




12/2012 


12/2012 



galaxy images (©. We provide in a description of 
the content of all published materials. 

2. HST/WFC3 IMAGING MOSAICS 
2.1. Fields and Filters 

A full desc ription of the CAND ELS o bserving program 
is giv en by iGrogin et al.l fj201 If) and iKoekemoer et all 
(|20ll . CANDELS is a WFC3 and parallel ACS, 902 
orbit HST imaging survey; here we concentrate on the 
WFC3 data only, which covers 800 arcmin 2 and is dis- 
tributed over five widely separated fields. There is a 
'deep' and a 'wide' component. About 125 arcmin 2 , di- 
vided over two fields, will have ~13 orbits per tile di- 
vided over three filters (F105W, F125W, and F160W), 
and the remaining area, distributed over five fields, will 
have 2 — 3 orbits per tile divided over two filters (F125W 
and F160W). These exposures reach -Hfi60W ~28 and 
Hfwow ~27 5ct magnitude limits for point sources in 
each filter for the 'deep' and 'wide' imaging, respectively. 

Along with this paper we electronically release struc- 
tural parameter catalogs for the three fields for which 
observations have been completed at the time of writing, 
as summarized in Table [H namely, th e Ultra Deep Sur- 
vey field fUDS. [Lawrence et al.|[2007[). the Cosmologica l 
Evolution Survey field (COSMOS lScoville et all [2007]) 
(both 9' x 24' arcmin and each at 'wide' depth), and 
the Great O bservatories Origins De ep Survey-South field 
fGOODS-S. rGiavalisco et"at|[200l : 'wide' over 4' x 10'; 
'deep' over 7' x 10'. Note that these dimensions refer to 
the CANDELS coverage, not to the dimensions of the 
data sets that define the original surveys. 

The CANDELS observations are augmented by previ- 
ously obtained WFC3/I R data from the Early Release 
Science program (ERS, IWindhorst et~aT1 MnS) in the 
Northern part of the GOODS-S field (4' x 9' at 2-orbit 
depth in F098M, F125W, a nd F160W) and the Ultra 
Deep Field (UDF) program (jBouwens et alj [2010) cm- 
bedded in the GOODS-S deep area (1 pointing with 
-15 orbits in F105W and F125W, and 28 orbits in 
F160W). The electronically available catalogs published 
along with this paper, derived from the finalized CAN- 
DELS imaging products at the date of publication, will 
be updated upon completion of the two remaining fields: 
GOODS-North a nd the Extended Groth Strip (EGS, 
iDavis et all 120071 ). Image, weight (inverse variance), 
and exposure time mosaics are prepared by drizzling 
the individual exposures onto a grid with rescaled pixel 
sizes of 0.06" . T h is pro cedure is described in detail by 
IKoekemoer et all (|2011f) . 

2.2. Object Identification 
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Figure 1. Left: Comparison between two HST/WFC3 F160W model PSFs. In red we show the TinyTim model PSF, in black a hybrid 
of a stacked star (outside a radius of 3 drizzled pixels, i.e. , 0.18") and the TinyTim model (inside that radius). The two sets of curves 
show the enclosed flux (left-hand y-axis) within the radius r shown on the x-axis, and the surface brightness per pixel (right-hand axis) at 
radius r. The vertical lines indicate the original WFC3 pixel scale (0.13") and the drizzled pixel scale of the mosaics used for the profile 
fitting in the paper (0.06"). The growth curves are normalized to the total flux of the model PSF within 10" and the growth curves of the 
empirical PSFs is forced to coincide with the that of the model PSF at 3". The surface brightness curves are azimuthally averaged and in 
units of fractional PSF flux per 0.06" pixel. Right: The residuals from individual stars after subtracting the TinyTim model PSF (red) and 
the hybrid model PSF (black). The systematic negative trend for the TinyTim model indicates that a larger fraction of the flux of a star 
resides at radii 0.2 — l"than what the TinyTim model predicts. Our hybrid model PSF does not show such a trend and follows the light 
profiles of the isolated stars out to radii of 3". 



We use Ei modif i ed version of 

SExtractor (jBertin fc Arnoutsl I1996T ) v2.5 to iden- 
tify objects in the CANDELS F160W mosaics. We refer 
to Galametz et al. (in prep.) for a full description of 
the source extraction but the main steps are as follows. 
The first modification to SExtractor is that we added 
a buffer between the isophotal area of each object and 
pixels used to estimate the local background. Note 
that this modification affects object identification and 
segmentation map construction, but the newly measured 
background is not used in the surface brightness profile 
fits (see below). The second modification is the removal 
of a bug that previously allowed disconnected regions 
to have the same value in the segmentation map, and 
therefore have the same identification number. Our 
modified version produces objects consisting of adjacent 
pixels only. 

SExtractor is run in dual mode, using the F160W 
mosaics both as the detection images and the measure- 
ment images. Indeed we found slight differences with re- 
spect to single-mode SExtractor output, which we wish 
to avoid in order to allow direct comparisons with dual- 
mode photometry of objects detected in F160W and mea- 
sured in other filters. 

Catalogs for two sets of SExtractor detection param- 
eters, the 'cold' setup and the 'hot' setup are produced 
sep arately and th e n com bined in an approach introduced 
by iBarden et al.l ((2012). The 'cold' setup focuses on 
optimal segmentation of relatively bright objects. This 
avoids that star-forming galaxies are spuriously split up 
into multiple components by deblending too aggressively. 
The 'cold' setup selects objects with 0.75a detections 
over 5 adjacent pixels after smoothing with a top-hat 
kernel with a diameter of 9 pixels. Objects arc de- 



blended adopting 9 logarithmic sub-thresholds relative 
to the maximum count value of an object and using a 
minimum contrast of 0.0001. Given that the PSF has a 
FWHM much smaller than that of the smoothing kernel 
used in the 'cold' setup, this approach is suboptimal for 
detecting small, faint objects. 

This is mitigated by the 'hot' setup, which uses a Gaus- 
sian smoothing kernel with a FWHM of 4 pixels that is 
similar to that of the PSF and selects objects with 10 ad- 
jacent pixels with 0.7a fluxes. Deblending is done with 
the default SExtractor parameters: a minimum contrast 
of 0.001 and 64 logarithmic sub-thresholds. 

Each object in the 'hot' catalog that falls within an 
aperture around an object in the 'cold' catalog is con- 
sidered part of the latter and removed from the com- 
bined catalog. This elliptical aperture is 2.5 times 
the elliptical Kron radius as measured with the default 
SExtractor parameters. The size and shape of the el- 
lipse are calculated from the second-order moment of the 
flux distribution. In the combined segmentation map the 
segmented pixels of the removed 'hot' object are added 
to the object in the 'cold' segmentation map. 

In total, we identify 109,533 unique objects in the 
F160W mosaics. HST resolution enables us to resolve 
close companions, which is essential when one is inter- 
ested in the measurement of the structures of individ- 
ual galaxies. 7% (9%) of objects in CANDELS 'wide' 
('deep') have one or more neighbors within 1"; they 
would typically not be resolved in ground-based imag- 
ing data sets. 

The SExtractor measurements of object magni- 
tudes, size, shape and orientation are used to provide 
GALFIT with initial guesses for the fitting parameters (see 
£ |4.2p . This helps in terms of computing time and ensures 
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that GALFIT converges to the global x 2 minimum in pa- 
rameter space. 

2.3. Background Estimation 

The image cutouts we use for model fitting ( 34.1[) arc 
too small for GALFIT to optimally determine the back- 
ground flux level. In order for the image to be large 
enough for an accurate, simmulataneous determination 
of both the structural parameters and the background 
flux level, the number of objects that would have to be fit 
simultaneously would become impractical, if not impos- 
sible in some cases (at present, GALFIT cannot fit more 
than 110 objects at one time). 

The pipeline GALAPAGOS (|Barden et alll2012f ) reme- 
dies this problem. O riginally written to analyze ACS 
mosaics from GEMS (jRix et all 120041) . GALAPAGOS de- 
termines the background from the full mosaic and then 
runs GALFIT using this background value along with 
initial guesses for the input parameters based on the 
SExtractor output described above. 

The independent background estimate is the main fea- 
ture of GALAPAGOS that elevates it beyond the level of 
'just' a smart wrapper for SExtractor and GALFIT . In 
short, ignoring pixels within a certain distance of an ob- 
ject, the flux is computed in a series of annuli, searching 
for a converging flux level at sufficiently large distances 
from the target object. This background flux level is pro- 
vided to GALFIT , which keeps it fixed at that value. In 
£15.11 we investigate the contribution of the uncertainty 
in the background to the uncertainties in the structural 
parameters. 

3. POINT SPREAD FUNCTION 

Besides an accurate noise model and background esti- 
mation, good knowledge of the PSF is essential for the 
accuracy of the parameter estimation. This necessitates 
the production of custom-made PSF models for each of 
the mosaics and each of the filters. These PSF mod- 
els should be used for all PSF-dependent analyses of the 
public CANDELS imaging mosaics, and are therefore re- 
leased here (see $6] for a summary of the published ma- 
terial) . 

The PSF has a FWHM of -0.17" for the F160W filter. 
Individual exposures (typically 800s) are dithered to im- 
prove on the originally coarse sampling - the pixel scale of 
the WFC3/IR ca mera is 0.13". Using the TinyTim pack- 
age (jKristl H995I ) we construct a model PSF for the re- 
spective imaging mosaics for the different fields and fil- 
ters. Briefly, TinyTim PSF models are created for the 
center of the WFC3 detector assuming a G2V spectral 
type. The PSFs are sub-sampled by a factor of ten 
in order to aid in aligning them with the CANDELS 
dither pattern. They are then re-sampled to the WFC3 
pixel scale and a kernel is applied to replicate the effect s 
of inter-pixel capacitance ([Hubert &: McCullougMl2l)lTI) . 
Finally, PSFs at each dither position are distortion- 
corrected and combined with the same drizzle parameters 
used in producing the imaging mosaics (|Koekemoer et all 

As an example we show the growth curve and surface 
brightness profile of the F160W PSF for the UDS (Figure 

13 http:/ /astro-staff.uibk.ac.at/~m.barden/galapagos/ 



Q]- red lines in the left-hand panel). To verify the accu- 
racy of this model PSF we subtract it from a sample of 46 
isolated stars in the F160W UDS mosaic (after normaliz- 
ing the star fluxes within a radius of 3", and examine the 
residual enclosed flux (Figure [1] - red lines in the right- 
hand panel). The systematic residual of — —4% seen 
at —0.2" implies that the model PSF contains —4% less 
light outside an aperture — 0.2" than point sources in the 
act ual d ata. Background levels, computed as described 
in §2.3[ are low compared to the stellar fluxes, even at 
radii of 1 — 3", such that the enclosed flux is affected by 
less than 1% at any radius. In other words, background 
errors cannot explain the large differences between stellar 
and model PSF growth curves at relatively small radii. 
Such an erroneous PSF model, which is seen for other 
fields and filters as well, can affect the measurement of 
structural parameters of small, concentrated objects at 
a level that exceeds the formal, random uncertainty. 

The difference between the stars and the 
TinyTim model occurs in relatively 'dark' areas 
seen in Figure [2] The center, the diffraction spikes, and 
the dot-like features associated with the first Airy ring 
are well reproduced by TinyTim , whereas regions in 
between those features are too dark. The origin of this 
feature remains unexplained. 

We use the 46 isolated stars shown in the right-hand 
panel of Figure [1] in median-stacked form as an alter- 
native PSF representation. The problem with this ap- 
proach is the variety in the sub-pixel positioning of the 
stellar images, which leads to broadening in the PSF 
model when stacking a number of stars. 

In order to provide an accurate PSF model at all radii 
we take the median-stacked star and replace the central 
pixels (within a radius of 3 pixels from the center) by 
the TinyTim model PSF. The flux values of these pix- 
els are normalized such that the total flux of the newly 
constructed hybrid PSF model is the same as that of the 
stacked star. The accuracy of the model PSF in the cen- 
tral region is confirmed by its good correspondence with 
the flux distributions of the small number stars that hap- 
pen to be well- aligned with the brightest pixel. Note that 
there are too few such stars to produce a stacked PSF 
with sufficient signal-to-noise ratio at larger radii. 

The growth curves, surface brightness profiles and 
residuals are shown in black in Figure [TJ The residu- 
als show much smaller systematic effects than the best- 
effort TinyTim model PSF, indicating that the hybrid 
PSF provides a good description of the growth curves of 
individual stars out to 3". 

We construct equivalent hybrid PSF models for all 
fields and filters. For the UDS and COSMOS fields we 
created PSF models for the F125W and F160W PSF fil- 
ters. For GOODS-S the situation is more complex. For 
the ERS region we have F098M, F125W and F160W PSF 
models; for the 'wide' and 'deep' regions we created sep- 
arate F105W, F125W and F160W PSF models; for the 
UDF, which is very small and contains few stars, we sim- 
ply make use of a single star. Replacing the central pixels 
by model values does not work well for a single star, but 
the chosen star is bright and isolated, and happens to 
be be relatively well-centered on the brightest pixel, and 
thus serves its purpose well. 

The fact that variations are < 2% for all stars implies 
that PSF differences due to variations in position, color, 
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Figure 2. Left: Ratio of our hybrid PSF and the TinyTim model, both for the UDS field, with dark regions lower levels in the TinyTim model 
(the darkest regions have values > 2). The image is 3.6" x3. 6", and the circle, which has been included to guide the eye, has a radius of 
0.6". The area near the center has a constant value by construction as we use the TinyTim model in our hybrid for the central pixels. We 
show our hybrid PSF with two different stretches in order to to emphasize the bright, central region (Middle) and the faint, outer region 
(Right). At each radius, the brightest parts of the PSF are reproduced quite accurately by the TinyTim model; rather, it is in between the 
diffraction spikes and Airy rings that the deviations are largest, and at a level that cannot be explained by noise. These features cause the 
systematic differences between the light profiles of stars and the TinyTim model shown in Figure [l] 



and magnitude are negligible for our GALFIT measure- 
ments. We note that small variations in dither pattern, 
orientation and general mosaic geometry exist within the 
distinct fields. This implies that some fraction of the 
objects will not have precisely the correct PSF model. 
This is unavoidable when using stars to produce the PSF 
model. 

4. MODELING LIGHT PROFILES 
4.1. Preparation of Image and Noise Cutouts 

For each of the 109,533 objects, an image cutout, in 
units of electrons per second, is produced from the large 
mosaic. The size of each rectangular cutout is deter- 
mined by the Kron radius, ellipticity and position angle 
as measured by SExtractor , and is chosen to enclose 
an ellipse with major axis length 2.5 times the Kron 
ellipse (Aim age x KRON radius m SExtractor par- 
lance). This typically corresponds to ~ 30 times the 
half-light radius of an object. 

A noise map of the same size is also produced. The 
pixel values in the noise map have units of electrons per 
second, and are the square root of the total variance. 
The total variance is estimate d by starting with the co m- 
puted variance, referred to by Kockemoer et al.l (|2011[ ) as 
the 'intrinsic' term (e.g., background noise and readout 
noise), in units of electrons. We add the variance at each 
pixel from the Poisson noise due to the objects them- 
selves. We divide by the (computed) exposure time for 
each pixel to make the units consistent with those of the 
images. 

The images and noise maps are used by GALFIT as 
inpu. We note that GALFIT does not have all necessary 
information about the image characteristics required for 
a self-contained noise-estimate as a result of the extensive 
process from raw, single exposures, to final, drizzled mo- 
saics. Therefore, especially for bright, compact objects 
it is essential to use the noise maps we construct rather 
than the GALFIT 's own noise estimation. In practice, for 
objects brighter than -£Tfi60W ^ 22 — 23 - correspond- 
ing to the typical background WFC3/F160W flux level 
measured over an area the size of a typical object - the 
extrinsic noise term begins to dominate. 



We compute signal-to-noise ratio (S/N) for each ob- 
ject from their images and noise maps, integrating over 
the pixels that belong to the objects according to the 
segmentation map. An object with i?Fi60W ~ 22 has 
S/N ~ 100 - 200 for the 'wide' data. These S/N es- 
timates are used in our derivation of the measurement 
uncertainties in |J5J 

4.2. GALFIT Setup 

The image cutouts, along with the noise model, the 
appropriate PSF model, and the pre-determined, fixed 
background level, are provided to GALFIT , which is used 
to find the best-fitting Sersic model for each object. The 
fitting parameters are total magnitude (m), half-light ra- 
dius (r) measured along the major axis, Sersic index (n), 
axis ratio (q), position angle (PA), and central position. 
Initial guesses for these parameters are taken from the 
SExtractor detection catalog. A constraints file is con- 
structed so that GALFIT is forced to keep the Sersic in- 
dex between 0.2 and 8, the effective radius between 0.3 
and 400 pixels, the axis ratio between 0.0001 and 1, the 
magnitude between and 40, and, between -3 and +3 
magnitudes from the input value (the SExtractor mag- 
nitude). 

Neighboring objects in each image cutout are fit simul- 
taneously or masked out, depending on their brightness 
compared to the main target: galaxies are fit simulta- 
neously if they are less than 4 magnitudes fainter than 
the main target; stars are fit simultaneously if they are 
less than 2 magnitudes fainter. In order to produce a 
good model, it is sometimes necessary to also fit ob- 
jects outside the image cutout. The decision to do so 
depends on the contribution of objects outside the image 
cutout to the flux in the cutout. SExtractor or previ- 
ously obtained GALFIT measurements are used to make 
informed decisions. The entire decision process, carried 
ou t by GALAPAG O S , is quite sophisticated and we refer 
to iBarden et al.l (|2012l ) for a full description. Note that 
the segmentation maps described in H2.2l are only used to 
identify objects, and not to identify the pixels that are 
used in the Sersic fits. 

For illustrative purposes, the images of 10 galaxies 
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Figure 3. The distribution of the GALFIT -derived parameters F160W magnitude, Sersic index, effective radius, axis ratio, and position 
angle for objects in the three CANDELS fields analyzed here (COSMOS, GOODS-S outh , and UDS). The panels with gray scale distributions 
represent the objects with flag value zero (these are objects with good fits - see i]4,3j l. On the diagonal we show for the flag zero sample 
the histogram of each of the five parameters in units as shown on the x-axis. The typical galaxy has magnitude -Hfi60W ~ 25, n ~ 1, 
r e ~ 0.3", and q ~ 0.4. The colored panels show the distribution of ob jects with flag value one (red) and flag value two (blue) — these are 
objects with, respectively, suspicious and bad fits as explained in £14.31 



and their fits are shown in Figure [4] All best-fitting 
GALFIT parameters are listed in Table [H and shown in 
Figure [3] The typical galaxy with a good Sersic model 
fit (see §4.31 for the definition of 'good') has magni- 
tude i?Fi60W ~ 25, Sersic index n ~ 1, effective radius 
r e ~ 0.3", and axis ratio q ~ 0.4. The position angle 
distribution is not entirely uniform; this is due to un- 
resolved sources, mostly stars. For resolved sources the 
distribution is uniform. 



4.3. Flags 

In Table [2] we flag objects with suspicious (flag value 
1), bad (2), or non-existent (3) fitting results. All other 
objects (good fits) have flag value of zero. Suspicious 
fitting results (with flag value 1) are those where the 
GALFIT magnitude deviates from the expected magni- 
tude by more than three times the uncertainty in the 
GALFIT magnitude derived as described below, in fcl5.ll 
This comparison is illustrated in Figure [5J where we 
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Table 2 

GALFIT fitting results and derived uncertainties. This table gives a sample of the UDS F160W catalog. Electronic versions of the full 
catalogs for all fields and (near-infrared) filters, as listed in Table [l] are published online (see Jj6] m is the AB magnitude, r is the 
half-light semi-major axis in arcseconds, n is the Sersic index, q is the projected axis ratio, PA is the position angle in degrees {PA = 
corresponds to North; PA = 90 corresponds to East). S/N is the signal-to-noise ratio measured over the pixels in the segme ntati on map 
that are attributed to the object. The listed uncertainties are based on population statistics and derived as described in i|5,ll and 
correspond to the half-width of the 68%-confidence interval (i.e., these are la error bars in the case of gaussianity) . 



ID 


RA 


DEC 


FLAG 


m ± <5m 


r ± 5r 


n ± 8n 


q±Sq 


PA ± SPA 


S/N 




J2000 


J2000 




AB mag 








dcg 




1 


34.223766 


-5.278053 


3 


-999 ± -999 


-999 ± -999 


-999 ± -999 


-999 ± -999 


-999 ± -999 


5.43 


2 


34.223904 


-5.277949 





24.41 ± 0.40 


0.43 ± 0.20 


0.80 ± 1.03 


0.18 ± 0.09 


41.61 ± 14.32 


6.88 


3 


34.223492 


-5.277952 





25.17 ± 0.24 


0.18 ± 0.06 


0.27 ± 0.28 


0.74 ± 0.37 


-52.04 ± 23.24 


8.57 


4 


34.265106 


-5.277749 


1 


23.06 ± 0.19 


0.27 ± 0.07 


1.49 ± 0.84 


0.72 ± 0.19 


67.92 ± 10.47 


13.33 


5 


34.295372 


-5.277648 


1 


21.06 ± 0.09 


1.16 ± 0.11 


1.55 ± 0.28 


0.80 ± 0.08 


-62.85 ± 3.57 


33.46 



Table 3 

Systematic and random uncertainties in structural parameters for CANDELS 'wi de' i maging in F160W. The systematic uncertainties are 

derived from GALAPAGOS processing of simulated images as described (se e !|5.2t ; the random uncertainties are inferred from 
GALAPAGOS processing of different data sets of the same galaxies as described in £|5.1I . The systematic uncertainties are given first, followed 
by random uncertainties as error bars. For m the units are in magnitudes; for the other parameters the uncertainties are relative, in linear 
units (they correspond to percentages). The galaxy samples are split by Sersic index (at n = 3) and effective radius (at r e = 0.3"). 







Am 






An 






m 


n < 3 


n > 3 


r c < 0.3" 


r e > 0.3" 


n < 3 


n > 3 


r e < 0.3" 


r e > 0.3" 


21 


-0.02 ± 0.01 


-0.01 ± 0.01 


-0.02 ± 0.01 


-0.01 ± 0.01 


-0.01 ± 0.02 


-0.03 ±0.03 


-0.02 ± 0.02 


-0.00 ± 0.02 


22 


-0.01 ± 0.02 


0.00 ± 0.02 


-0.01 ± 0.01 


-0.01 ± 0.02 


-0.03 ± 0.04 


-0.06 ± 0.06 


-0.05 ± 0.04 


-0.01 ± 0.05 


23 


-0.01 ± 0.03 


0.06 ± 0.07 


-0.01 ± 0.03 


0.01 ± 0.05 


-0.01 ± 0.10 


-0.09 ± 0.16 


-0.01 ± 0.09 


-0.03 ± 0.13 


24 


-0.01 ± 0.06 


0.13 ± 0.13 


-0.01 ± 0.05 


0.05 ± 0.12 


-0.01 ± 0.24 


-0.25 ± 0.33 


-0.01 ± 0.23 


-0.11 ± 0.33 


25 


-0.04 ± 0.13 


0.10 ± 0.21 


-0.03 ± 0.12 


0.03 ± 0.25 


0.02 ± 0.49 


-0.49 ± 0.58 


0.00 ± 0.45 


-0.20 ± 0.68 


26 


-0.14 ± 0.26 


0.02 ± 0.31 


-0.13 ± 0.25 


-0.09 ± 0.44 


0.31 ± 0.86 


-0.84 ± 0.92 


0.23 ± 0.84 


-0.16 ± 1.21 


27 


-0.47 ±0.41 


-0.37 ± 0.54 


-0.46 ± 0.42 




0.99 ± 1.28 


-1.59 ± 1.46 


0.87 ± 1.29 








Ar 








Aq 






m 


n < 3 


n > 3 


r c < 0.3" 


r c > 0.3" 


n < 3 


n > 3 


r c < 0.3" 


r c > 0.3" 


21 


0.00 ± 0.01 


0.00 ± 0.016 


0.01 ± 0.01 


0.00 ± 0.01 


0.01 ± 0.01 


-0.01 ± 0.02 


0.01 ± 0.01 


0.0 ± 0.01 


22 


0.01 ± 0.02 


-0.02 ± 0.04 


0.02 ± 0.02 


-0.01 ± 0.02 


0.00 ± 0.02 


0.01 ± 0.03 


0.00 ± 0.02 


0.0 ± 0.03 


23 


0.00 ± 0.04 


-0.10 ± 0.11 


0.00 ± 0.03 


-0.03 ± 0.06 


0.00 ± 0.05 


-0.01 ± 0.08 


0.00 ± 0.04 


0.00 ± 0.06 


24 


0.01 ± 0.08 


-0.22 ± 0.19 


0.01 ± 0.08 


-0.10 ± 0.15 


-0.02 ± 0.11 


-0.03 ± 0.18 


-0.02 ± 0.10 


-0.02 ± 0.15 


25 


0.04 ± 0.19 


-0.19 ± 0.33 


0.04 ± 0.18 


-0.09 ± 0.33 


-0.07 ± 0.25 


-0.07 ± 0.35 


-0.07 ± 0.24 


-0.04 ± 0.35 


26 


0.12 ± 0.43 


-0.11 ± 0.55 


0.12 ± 0.42 


-0.11 ± 0.63 


-0.17 ± 0.51 


-0.09 ± 0.59 


-0.17± 0.51 


0.06 ± 0.69 


27 


0.27 ± 0.75 


0.22 ± 0.85 


0.27 ± 0.76 




-0.28 ± 0.81 


0.40 ± 0.90 


-0.25 ± 0.82 





show the GALFIT and SExtractor magnitudes for the 
UDS F160W measurements. The expected magnitude 
is the BEST magnitude from SExtractor , corrected 
for the systematic, magnitude-dependent offset between 
GALFIT model magnitude and the BEST magnitude, in- 
dicated by the red line in the figure. The BEST mag- 
nitude is measured in i?Fi60W, adding the color, e.g., 
-/fi25W — ^Fi60Wj measured over the £?fi60W segmen- 
tation map. This offset between the SExtractor and 
GALFIT magnitudes is 0.1 mag for bright sources and 
increases to 0.3 mag for faint sources (see Figure [5]). 
This offset has bee n noted by many auth ors before (e.g., 
Holden et afl 120051: iBlakeslee et al 1120061 : lHaussler et al I 
20071) . The distribution of GALFIT fitting parameters 
with flag value I is shown in red in Figure [Sj 

Objects with flag value 1 are not necessarily bad fits 
and in some cases can be used without problem. How- 
ever, we recommend to assess on an object-by-object ba- 
sis whether the results can be used and to examine in 
those cases the GALFIT model and residual, which are 
also released along with this paper (©. 

Bad fitting results (with flag value 2) are those for 
which one or more parameter reached the constraint 
value forced onto GALFIT . Note that these fits may in 



fact truly represent the best possible Sersic profile, but 
that the inferred structural parameters are in most cases 
not astrophysically meaningful. This should be assessed 
on an object-by-object basis and refined, hand-tuned fit- 
ting. The user may adopt inferred structural parameters 
at their own risk. Objects that are fit simultaneously 
with the target galaxy are allowed to reach those con- 
straint values - the flag value for the target remains 
in this case. A flag value 2 is also assigned in case the 
axis ratio drops below 0.1, which GALFIT tends to do in 
the case of seemingly converged, yet bad fits of generally 
faint objects. The distribution of GALFIT fitting parame- 
ters of objects with flag value is shown in blue in Figure 
[3J Finally, non-existent results (with flag value 3) are 
simply those fits that 'bombed', in GALFIT parlance. 

The depth of our catalog of structural parameters for 
objects with flag value zero is i?Fi60W ~ 25 as can be 
seen in the -Hfi60W magnitude histogram in Figure [H 
Accuracy and precision of the measurements for flag-zero 
objects are discussed at length in £J31 but this canonical 
depth of i?Fi60W ~ 25 is two magnitu des shallower than 
the 5ct limit for CANDELS 'wide' (jKoekemoer et al.l 
I2011h . Thus, for objects in the range 25 < i?Fi60W < 27 
that are securely detected by CANDELS 'wide' a mea- 
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Figure 4. Ten galaxy images (left) from the UDS F160W mo- 
saic and their GALFIT models (middle) and residuals (right). These 
objects have been chosen to have a range in magnitude (increasing 
from -ffpigow ~ 20 at the top left to -Hfi60W ~ 25 at the bottom 
right) and structure. UDS-7335 and UDS-2813 have n ~ 4 and 
smooth light profiles characteristic of early- type galaxies; UDS- 
24786 and UDS-7373 have n ~ 1 and show more substructure, 
characteristic of star-forming, late-type galaxies. For the faintest 
galaxies shown here the GALFIT model still captures the basic struc- 
tural parameters such as size and shape. In several cases multiple 
obj ects are fit simultaneously, illustrating the procedure described 
in fl4~2l 

surement of even their basic structural parameters will 
typically not be possible. There is a population of flag- 
zero objects that extends to fainter magnitudes, up to 
-Hfi60W ~ 28 (most easily seen in the magnitude-radius 
panel of Figure [3]); these objects are located in the UDF. 



5. MEASUREMENT UNCERTAINTIES 



5.1. 



Random Uncertainty Estimates from Internal 
Comparison 

A robust way to assign random uncertainties to the 
GALFIT measurements is to compare measurements of 
the same objects in different data sets. For the ob- 
jects detected in the 'deep' region of GOODS-S, we re- 
run GALAPAGOS on a shallower mosaic that has the same 
depth as the 'wide' CANDELS imaging. All steps are 
performed in a manner that is precisely analogous to 
that applied to other data sets. First, the appropriate 
hybrid PSF model is constructed as described in SjSj and 
GALAPAGOS is used to estimate the background level, as 
described in 12.31 before running GALFIT . The input seg- 
mentation map and SExtractor catalogs are identical 
for the for the 'deep' and the 'shallow' versions of the 
mosaics. Therefore, neighboring objects are treated in 
exactly the same way in both cases, and differences in 
GALFIT fitting results are entirely driven by the noise 
and background properties of the images such that the 
variation between the two sets of measurements reflect 
the uncertainty. Note that when running GALFIT we do 
not distinguish between pixels included and excluded by 
the segmentation map, which is only used to identify ob- 
jects. 

The measurement uncertainties in the structural pa- 
rameters depend foremost on the S/N, as can be seen 




Figure 5. Comparison of GALFIT and SExtractor (BEST) F160W 
magnitudes for objects in the UDS. The y-axis shows the difference 
between the two. Black points represent objects that we have as- 
signed flag value (good fits); red points have flag value 1 (suspect 
fits); magenta points have flag value 2 (bad fit s - b right stars have 
bad fits because they are point sources). See i]4.3l for a full expla- 
nation of the flag definitions. The red line is a running median 
for the black and red points, quantifying the systematic offset be- 
tween GALFIT and SExtractor total magnitudes. Note that these 
systematic offsets are in excess of the random uncertainty in the 
GALFIT -based total magnitudes. 



0.5 - 



-0.5 
0.2 



oo o 
o 



-0.2 - 



0.2 




1000 



S/N 



Figure 6. Difference between the structural parameter estimates 
from deep and shallow images of 6492 objects in GOODS-S as a 
function of S/N in the shallow up. The sprea d re flects the uncer- 
tainty in the measurements, which we use in £|5.1l to assign uncer- 
tainties to the structural parameter measurements of all galaxies. 
Although dependencies on other parameters exist (see Figures [7] 
and[$), the uncertainty depends to first order on S/N. 
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Figu re 7. Cor relations among structural parameter measurement uncertainties (also see, e.g., IHiiussler et~aTl l2007; Guo eTaLl |2009| ; 
Bruce ct al. 2012). All quantities A arc the difference between the structural parameters are measured from the deep and shallow images of 
6492 galaxies. A is normalized to S/N = 50 as described in the text; that is, first-order effects on S/N are removed. Particularly relevant 
are the strong correlations between the uncertainties in magnitude (m), Sersic index (n) and effective radius (r): uncertainties in m and 
r correlate such that any overestimate (underestimate) in r corresponds linearly with an overestimate (underestimate) in the total flux; 
analogously, if the total flux is overestimated (underestimated) by 5%, the Sersic index is typically overestimated (underestimated) by as 
much as 25%. The shape and amplitude of these 'error ellipses' (or rather, 'ellipsoids', as they are at least 3-dimcnsional) depend on the 
parameters themselves and therefore do not represent the covariance matrix of the uncertainties for the sample as a whole. 
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Figure 8. Correlation between structural parameter measurement uncertainties and measured structural parameters. On the x-axis 
three parameters (magnitude, m, Sersic index, n, and effective radius, r) as measured from the GSD deep region are shown. The galaxies 
displayed in the first column of panels (showing m), are chosen within narrow ranges of n and r around the median values in the full catalog. 
In the second (third) column the same is done for n (r): m and r (n) are chosen within a narrow range. On the y-axis the difference A 
between the measurements from the deep and shallow images of the same objects are shown: the scatter in A reflects the uncertainty, 
indicated in red by the running 16-84%-tile ranges. As in Figure[7] all A are normalized to S/N = 50. Strong second-order effects remain, 
most clearly seen in the dependence of the uncertainty in r and m on n and the dependence of the uncertainty in n and m on r. 



in Figure |6l However, there are several complicating fac- 
tors that significantly affect the true uncertainties: 1) the 
uncertainties in different parameters are correlated (see 
Figure [7]); 2) after normalizing for S/N, the uncertain- 
ties themselves also depend on the parameter values in 
a non-trivial manner (see Figure [8] for a demonstration), 
and; 3) the covariance between the measurement uncer- 
tainties changes with the values of the parameters and 
their respective uncertainties. 

We adopt an empirical approach in which the two fit- 
ting results (from the deep and the shallow data) of sets 
of similar galaxies are used to calculate, throughout pa- 
rameter space, the uncertainties in each of the measure- 
ment parameters. We assume that the measurement un- 
certainties depend on m, n, and r (see Figurc|5]), but not 
on any other parameter: we do not see evidence for a 
correlation between q (or, more obviously, PA) and the 



uncertainties in any of the parameters. 

We construct a parent sample of 6492 galaxies g with 
GALFIT measurements from deep and shallow data and 
a flag value less than two (see Sec I4.3[) for both fits. For 
each galaxy gi we identify the 200 most similar galaxies 
in the parent sample g by computing the normalized dis- 
tances pi — pj in the three-dimensional parameter space 
spanned by m, n, and r. That is, pi is defined as 

Pi = (mi/a(m), log raj a {log ft), log r { /a (log r)), (1) 

where a denotes the standard deviation in the respec- 
tive parameters. These factors are introduced in order 
to make differences in parameter values comparable and 
produce dimensionless, normalized distances. 
For the 200 galaxieiQ most similar to the target galaxy 

14 This number can be changed by a factor of two without sig- 
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gi, as defined by pi, we compute the 16-84%-tile range in 
the differences between the parameters inferred from the 
deep and the shallow imaging. This provides an uncer- 
tainty estimate for each parameter as a function of the 
parameters m, n, and r: 

Si = S{pl) = {Sm i: dlogm, Slogn, Sqi, SPAi) (2) 

Si is normalized through multiplying by S/Ni, as mea- 
sured from the segmentation map (see @; this removes 
the first-order dependence on S/N and allows us to com- 
pute the error budget for galaxies from images with dif- 
ferent exposure timeJ^l By repeating this computation 
for all galaxies g we map out the measurement uncer- 
tainties throughout parameter space sampled by p. 

This large matrix serves as a database which we use 
to estimate the uncertainties in the GALFIT parameters 
for all galaxies in all our fields. As before, we calcu- 
late the distances pi — Pj, where pj represents the ob- 
jects in the database, and pi represents the galaxies to 
which we want to assign measurement uncertainties. We 
take the average 5j of the 25 'nearest' galaxies gj in the 

databasePl and divide by S/Ni to provide Si with the ap- 
propriate amplitude. This quantity represents the mea- 
surement uncertainties in all parameters (mi, Uj, r^, 
and PAi). Note that the uncertainties are correlated 
with each other, approximately as shown in Figure [7l 
This figure serves as a mere illustration; not only the 
amplitude of the uncertainty but also the covariance de- 
pends on the measurement parameters themselves (i.e., 
on pi). 

All measurements and their uncertainties (marginal- 
ized over all other parameters are described above) are 
provided in Table [2j For uncertainties that are computed 
in logarithmic units we give uncertainties in the custom- 
ary linear units in the table, where Sn = ln(10) n <51ogn 
in the case of, for example, the Sersic index. The mag- 
nitude and its uncertainty are kept in the usual archaic 
form. For reference, average random uncertainties (as 
well as the systematic uncertainties - see below, in £j 5 . 2 1) 
are given in Table [3] for galaxies with different mag- 
nitudes, sizes, and Sersic indices. The bottom line is 
that for the CANDELS 'wide' imaging in F160W the 
parameters m, r e and q can be inferred with a ran- 
dom accuracy of 20% or better for galaxies brighter 
than i?Fi60W ~ 24.5, whereas n can be measured at 
the same level of accuracy for galaxies brighter than 
#Fi60W ~ 23.5. 

For those readers who wish to derive their own uncer- 
tainty estimates (for example, in a difference confidence 
interval, or for investigating asymmetric uncertainties) 
we provide our fitting results for the 'wide' sample in the 
GOODS-South, which consists of 6492 objects (Table HJ). 
The IDs of these objects match those in the F160W cat- 
alog for GOODS-South. 

While the uncertainty matrix is generated from F160W 
data, we can also use it to assign uncertainty estimates 
for the structural parameters measured in the other fil- 

nificantly changing the results. 

15 Here, we add the noise from the deep and shallow images in 
quadrature. Thus, rather than adopting the measurement from 
the deep image as 'truth', we compute the combined error from 
the two measurements. 

16 see Footnote Q3] 



Table 4 

GALFIT fitting results for the F160W 'wide' sample in 
GOODS-South (6492 objects) used for estimating the uncertainty 
esti mates through comparison with the 'deep' fitting results (see 
i|5.1l l. The ID matches the ID in the general catalog for 
GOODS-South. An electronic version of the full catalog is 
published online. For further explanation of the table entries, see 
Table |2] 



ID 


m 


r 


n 


1 


PA 


S/N 




AB mag 








deg 




2989 


24.09 


0.18 


1.16 


0.53 


29.84 


23.15 


3034 


24.05 


0.25 


0.56 


0.47 


-51.53 


23.71 


3113 


22.79 


2.12 


7.15 


0.53 


82.94 


20.68 


3148 


23.51 


0.16 


1.85 


0.27 


-68.60 


42.53 


3166 


25.17 


0.43 


0.71 


0.77 


24.50 


4.68 



ters, again using the S/N to provide the uncertainties 
with the correct amplitude. Several factors may com- 
promise this approach: the generally less smooth light 
distribution at shorter wavelengths will tend to increase 
in the uncertainties associated with single-component 
Sersic fits, while the smaller PSF will tend to decrease 
the uncertainty, especially for compact objects. How- 
ever, uncertainties depend to first order on S/N alone 
and there is no reason to assume that the correlations 
shown in Figures [7] and [8] substantially change with filter 
choice. 

We note that the reported uncertainties represent the 
convolved contributions from noise an d the unce r tainty 
in t he backgroun d leveh and we refer to lGuo et al.l (|2009f ) 
and iBruce et al.l ([2012) for detailed analyses of errors 
due to uncertainties in the background flux. To pro- 
vide an indication of the contribution of the uncertainty 
in the background we rerun GALFIT for 1000 randomly 
picked galaxies with estimated background values bi in 
the UDS, assigning background levels bj as calculated by 
GALAPAGOS for another 1000 randomly picked galaxies. 
The standard deviation in bj—bi is 1.4xl0~ 4 e/s/pix. We 
then rerun GALFIT to obtain 1000 sets of parameter esti- 
mates that we compare with the original estimates. Be- 
cause the variation in background is likely partially real, 
the variation in the structural parameter estimates pro- 
vides an upper limit on the contribution of background 
uncertainties to the measurement errors. This exercise 
shows that at most ~ 25 — 30% of the total error bud- 
get as derived above is due to uncertainties in the back- 
ground flux estimates. Hence, for most objects the accu- 
racy in the structural parameter estimates is limited by 
the S/N, and not the background flux level, even though 
the latter is not entirely negligible. Only for very faint 
sources (-Hfi60W > 25.5) with large sizes (r > 0.4") the 
uncertainty in the magnitude and structural parameter 
estimates starts to be dominated by the uncertainty in 
the background estimate, but for such faint sources the 
measurements are highly uncertain anyway (see £|5.2|) . 

5.2. Systematic Uncertainties from Simulated Images 

From the data itself it is difficult to infer system- 
atic uncertainties as the 'true' light distributions of the 
galaxies are unknown. We use simulated images with 
galaxies with known light distributions to quantify these. 
The simulated imag es are generated as described by 
lHaussler et al.1 (|2007f ). with the background noise taken 



12 



Structural Parameters of Galaxies in CANDELS 



directly from empty parts of the UDS F160W mosaic. 
The input Sersic profiles of the fake galaxies have the 
same magnitude, size and shape distributions as in the 
real data. The input profiles are convolved with the same 
PSF as used in the data analysis - this exercise is not de- 
signed to test for systematic effects due to errors in our 
PSF model. The simulated mosaic is then processed as 
described in Sj2] and the inferred structural parameters 
can be compared with the input values. 

In Table [3] we give the average difference between the 
input and output values of the various parameters, along 
with the sample average random uncertainties assigned 
as described above in N5.ll For the largest part of pa- 
rameter space, and especially for the magnitude range 
of interest for morphological studies of massive, high- 
rcdshift galaxies, the systematic differences are substan- 
tially smaller than the random uncertainties; that is, ran- 
dom uncertainties dominate. However, systematics are 
not negligible, especially for objects with large Sersic in- 
dex. Generally speaking, structural parameters can be 
measured with a precision and accuracy better than 15% 
down to -ffFi60W ~ 23. At fainter magnitudes, there 
are larger systematic and random uncertainties, espe- 
cially for large, high Sersic index objects. Since the typ- 
ical faint, high-redshift galaxies is small and has a low 
Sersic index, 10%-level accuracy and precision in the ba- 
sic size and shape parameters should be reached down to 
-^Fieow ~ 24.5. 

Although the tabulated values strictly apply to the 
'wide' imaging, the small changes per magnitude bin im- 
ply that the systematic effects to the 'deep' imaging are 
not substantially different. We encourage all users of the 
public catalogs to take the results given in Table [3] into 
account in their analysis. We do not correct the measure- 
ments for these systematic effects, as we cannot quantify 
the precise amount for each individual galaxy. The re- 
sults presented here serve as indications of the magnitude 
of the systematic uncertainties. Note that these system- 
atic uncertainties do not include errors in the WFC3 zero 
points or other uncertainties in the calibration. 

The variance in the difference between the input and 
output structural parameters reflect random uncertain- 
ties that are very similar to those inferred in H5.ll with 
a similar dependence on S/N and correlation between 
the various parameters. The reason that we choose to 
use real rather than simulated data to infer our random 
uncertainties is that single Sersic profiles do not (neces- 
sarily) represent the true light distribution of a galaxy 
very well, which could lead to underestimated uncertain- 
ties in the case that idealized simulated galaxies are used 
in the analysis. 

As mentioned above, the input structural parameters 
for the simulations are based on the observed distribu- 
tion of parameters. This introduces the risk that re- 
gions of parameter space with large systematic effects 
arc unjustifiably ignored by design. We test the de- 
gree to which very small galaxies can be recovered by 
our method by generating an additional set simulated of 
galaxy images. The input magnitudes are in the range 
23.5 < i?Fi60W < 24.5, corresponding to the faintest 
galaxies that we still have fairly precise measurements 
for (see Table [3]) . The structural parameters are chosen 
to lie on a linearly-spaced grid in (log(r), n, q, PA) 
space, with the semi-major axis r ranging from 0.3 to 30 



pixels, probing a much wider range than observed. 

We find that 95% of all simulated objects with semi- 
minor axis lengths of 0.3 pixels are recovered correctly, 
meaning that GALFIT does not crash or converges to an 
incorrect value, and that there are no systematic differ- 
ences between the input and recovered parameter values. 
We conclude that the lack of large numbers of such small 
galaxies in the real catalogs (as seen in, e.g., Figure [8]) is 
not due to limitations in spatial resolution. 

6. PUBLIC DATA RELEASE 

This paper is accompanied by the public release^ of 
a number of materials. Most importantly, we present 
catalogs containing 109,533 unique objects in three fields 
(GOODS-South, UDS, and COSMOS), with structural 
parameter measurements for each object in at least two 
filters, F160W and F125W, and, in some cases a third 
(F098m or F105W; see 32). Catalogs for two additional 
fields (GOODS-North and EGS) will be released upon 
completion of the survey, as indicated in Table [TJ 

The catalogs contain an identification number, coor- 
dinates measured from the F160W mosaics (i.e., the co- 
ordinates are the same in the catalogs for the different 
filters), a flag indicating the availability /reliability of the 
GALFIT model (JOJ), the structural parameters measure- 
ments (§3} and their random uncertainties ( §5.111 . A sam- 
ple of the UDS F160W catalog is given in Table H and 
the parameter distribution is shown in Figure |3] Image 
cutouts, model surface brightness profiles, and residuals 
from the model fit are also made available online, an ar- 
bitrarily chosen sub-sample of which are shown in Figure 
[4j F160W object segmentation maps (see §2.2p are made 
available for each field, and PSF models, constructed as 
described in §3j are published for each field and filter. 
These models combine stacked stars and TinyTim mod- 
els to provide optimal sampling and fidelity. 

The structural parameter estimates are not corrected 
for systematic uncertainties, which we quantify in §5.21 
and Table [3l but we encourage users to verify the impact 
of systematic errors on their analysis. Table [3] also con- 
tains average random uncertainties for a range of magni- 
tudes. This table serves as a guideline when choosing a 
magnitude limit, or verifying the accuracy and precision 
of the fitting results for a particular set of objects. 

This work is supported by HST grant GO-12060. Sup- 
port for Program number GO-12060 was provided by 
NASA through a grant from the Space Telescope Science 
Institute, which is operated by the Association of Uni- 
versities for Research in Astronomy, Incorporated, under 
NASA contract NAS5-26555. 
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